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'sj- ' pair production and explore analytically and numerically the solution of the 

■ non-Markovian kinetic equation. In the Markovian limit of the kinetic equa- 

QQ , tion we find an analytical solution for the single particle distribution function 

I of bosons and fermions. The numerical investigation for a homogeneous, con- 

I stant electric field shows an enhancement (bosons) or a suppression (fermions) 

of the pair creation rate according to the symmetry character of the produced 
particles. For strong fields non-Markovian effects are important while they 
. disappear for weak fields. Hence it is sufficient to apply the low density limit 

for weak fields but necessary to take into account memory effects for strong 
^ : fields. 



Abstract 



We analyze a quantum kinetic equation describing both boson and fermion 



^ . PACS: 12.38.Mh, 05.60.+W, 25.75.Dw 

I. INTRODUCTION 

A proper description of the pre-equilibrium evolution of the quark-gluon plasma [QGP], 
believed to be created in an ultrarelativistic heavy-ion collision can start from a transport 
equation that incorporates a source term and a collision integral [![]. The source term 
describes the production of pairs of particles and antiparticles while the collisions lead to 
thermalization. 

The formation of the QGP is assumed to proceed via the creation of a strong chromo- 
electric field in the region between the two receding nuclei after the collision. The field 
subsequently decays by emitting quark-antiquark pairs according to the nonperturbative 
tunneling process of the Schwinger mechanism 0-^]. Within the flux tube model a lot of 
promising research has been carried out . 
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The process of pair production within the Schwinger mechanism has been addressed by 
many authors P,|6|,p!0,15| in recent years, with the back reaction scenario also considered [16 



A recent apphcation of the Schwinger mechanism of pair creation to the QCD case, solving 
a transport equation in boost invariant variables with a simple collision and source term, 
has been performed in |T^. Another interesting calculation was provided by where the 
polarization of the created quark-antiquark pair in a string fragmentation model for inclusive 
reactions was discussed. 

While these studies have been very useful in exploring the preequilibrium physics of a 
heavy ion collision, there is still the open question how to link the field theoretical treatments 
with a kinetic theory. More recent investigations performed for QED have shown that a 
consistent field theoretical approach leads to a kinetic equation with a modified source term 
providing a non-Markovian evolution of the distribution function. This result was first 
obtained by Rau [0 for the case of fermion pair creation for a constant electric field using 
a projection method [^. Therein it was emphasized that the particle pairs are produced 
with a non-trivial momentum dependence and not only at zero longitudinal momentum as 
assumed in previous studies. 

A generalized treatment allowing for a time dependent field was given in pT|j2^ . Therein 
we derived a kinetic equation containing both boson and fermion pair production. For the 
case of a constant field our results agree with those of |lT9| , p3| . The properties of the source 
term itself such as the momentum dependence and the time structure have been studied 
2^^4[ for a constant and a time dependent field. A detailed analysis of the boson pair 



creation combined with a systematic numerical study of the time structure of the solution 
was provided in the approach of [^]. These studies are performed for weak fields where the 
Markovian approximation is valid. 

However, the appearence of non-Markovian aspects of the pair creation process which 
are expected for strong fields are not yet fully understood. Therefore the main goal of this 
article is the study of the time evolution of the system in a regime where memory effects 
become important. As an illustrative example we consider pair creation in QED. Attention 
will be paid to the numerical analysis of the evolution of the distribution function for both 
bosons and fermions. We explore the influence of the different quantum statistical properties 
of boson and fermion pair creation and how the differences depend on the strength of the 
external field. 

In Section II. A we discuss the kinetic equation with a new source term for pair production. 
In Section II. B we explore the numerical solution of the non-Markovian kinetic equation for 
weak and strong fields. We also study the influence of the symmetry character of fermions 
and bosons on the evolution of the distribution function. In Section II. C we discuss the 
analytic solution in the Markovian and the low density limits, respectively. We compare 
the results with the non-Markovian solution and discuss the differences. The results are 
summarized in Section III. Within this study we neglect any influence due to collisions as 
well as back reactions of the produced charged particles on the initial electric field. 



II. THE NON-MARKOVIAN KINETIC EQUATION 



2 



A. A source term with non-Markovian character 



We consider particle production in a strong external electric field which leads to an 
unstable vacuum that can decay by creation of electron-positron pairs. Using the field- 
theoretical model of charged particles in an external, homogeneous, time-dependent field 
characterized by the vector potential = (0, 0, 0, A(t)) with A(t) = A^it) and the resulting 
electric field E(t) = E^it) = —A(t) = —dA(t)/dt, the kinetic equation is derived starting 
from the Dirac- (Klein-Gordon-) equation for fermions (bosons). The transition from the 
in-state to the instantaneous, quasiparticle state at the time t has been achieved by a time- 
dependent Bogoliubov transformation. Details of this derivation are given in [^,^. As 
the final result we obtain the kinetic equation for the single particle distribution function 
f{P,t) =< 0\a'^p{t)ap{t)\0 > defined as the vacuum expectation value in the time dependent 
basis of creation and annihilation operators a^(t), ap{t) for electron states at the time t and 
the 3- momentum P 

^W±(t) /* rft'W±(tO[l±2/±(P,t')]cos[x(t',t)], 

where the upper sign (lower sign) in Eq. (|l|) corresponds to boson (fermion) pair cre- 
ation. Details of the derivation are given in [pT|-p3|. The momentum is defined as 
P = {pi,p2, P\\{t)), with the longitudinal momentum P\\it) = py — eA{t) where p\\ = p^. 
For fermion creation we find in agreement with 



(2, 



and for boson production our result agrees with 



w,,o ^ £^ ^ aw w_,o. (3) 



We define the total energy uj{t) = + P'^{t), the transverse energy = ^Jnl? and 
the transverse momentum p^ = {pi,P2)- Furthermore x{t',t) = 2[Q{t) — 6(t')] denotes the 
difference of the dynamical phases which are defined as 

e(t) = dt'uj{t') . (4) 



Equation (||) is characterized by the following properties, (i) The particles are produced not 
only at rest p\\ = as assumed in more phenomenological approaches, e.g. [|1^. In Fig. 1 the 
dependence of the source term in low density limit for fermion production on the parallel 
momentum 

Sl{P,t) = ^>V±(t) f rft'W±(t')cos[x(t',t)] , (5) 
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is shown. This calculation has been performed for a strong and weak constant electric 
field, E{t) = const. The production rate is peaked at about zero momentum, for positive 
momenta it approaches zero and for negative momenta it is dominated by oscillations due 
to the choice of a constant electric field. It increases with increasing strength of the external 



field. Similar results have been obtained recently by different authors [p!9|,pl|-p3| . (ii) The 
source term and the distribution function have a momentum dependence accounted for by 
the transverse energy e± and by the kinetic momentum P\\{t), i.e. once a solution for py = 
and = is obtained, the solution for nonvanishing momenta can be generated by simple 
variable transformations. Therefore, we drop the explicit notation of the dependence on P 
in the distribution function and in the source term. 

(iii) Furthermore, the kinetic equation (|l|) has non-Markovian character due to the ex- 
plicit dependence of the source term on the time evolution of the distribution function. The 
source term contains a time integration over the statistical factor [1 ± 2f±{P, t)] which may 
cause memory effects. This important property will be discussed within this article. Inves- 
tigating the differences of boson and fermion pair creation for weak and strong fields, we go 
beyond recent studies, e.g. by p2|,|23|, in the following subsections. 



B. Solution of the non-Markovian kinetic equation 

In the previous subsection we have discussed general features of the source term, now 
we want to study the time structure of the non-Markovian solution in detail by solving the 
kinetic equation 



df±it) 1 
dt 2 



W±{t) f rft'W±(t')[l±2/±(t')]cos[x(t',t)] (6) 

^— oo 



numerically. The non-Markovian character is obvious. The appearance of the statistical 
factor under the time integral means that the solution of the differential equation depends 
on the full time evolution of the distribution function and hence memory effects are included. 
This complicated structure requires a self-consistent scheme for the numerical solution. 
As initial conditions we use lim f±{t) = 0. We solve the differential equation within the 

t — >—oo 

standard methods and obtain the self-consistent solution by iteration. Another possibility to 
solve this equation is the direct solution by integration. Numerical investigations have shown 
that this method leads to the same result but more iteration steps are needed to find the 
solution and thus we decided to solve in its differential form. In order to demonstrate the 
numerical solution of the kinetic equation we choose the simple case of a constant electric 
field E{t) = E = const.. On the one hand this is the most simple possible case also 
for the numerical treatment since this assumption leads to a reduction of the number of 
integrations. On the other hand this is the standard Ansatz for the case that back reactions 
are not included |[T9|j23|| and therefore permits to compare the results. The vector potential 



is in this case 

AiT) = AiT)/e^ = -TE,/e, (7) 

where the dimensionless variable Eq = eE/e\ does not depend on time and the energy is 
given as 
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u,{r) = ^l + El{r-^TY. (8) 

In our calculations we keep the transverse momentum fixed and hence normalize in units of 
a constant e^. In Eqs. (0) and (H) we have introduced the dimensionless time variables 

r = t£±, (9) 

and 

Ar = Ar(p,.p,) = Mi|il. (10) 

This notation is also convenient to distinguish the weak field {Eq < 1) and strong field 
{Eq > 1) limits. 

For such a simple model case, we can solve the integral of the dynamical phase @ and 
obtain 



0(r) = ^(r-ArVo(r) + ^ln 



Eo{t - At) + cuoir] 



(11) 



In the following we discuss the numerical solution of the kinetic equation for bosons and 
fermions for different field strengths. In Fig. 2 we plot the distribution function for bosons for 
weak fields and strong fields. The particles are produced at about zero kinetic momentum. 
Because of the fact that we have no damping mechanism the distribution oscillates around 
a constant value. The oscillations to be seen in the figure are due to the choice of a constant 
electric field. In response to strong fields the frequencies of these oscillations increase while 
the amplitudes decrease compared with the limit at large times. Other Ansdtze for the time 
dependence of the external field such as a Gaussian shape |]22| may lead to a damping of 



these oscillations. But it is important to note that in order to describe a more realistic 
situation the back reaction of the produced particles on the initial field should be included 
16| , |23[| . The curves are normalized to the large time limit 



fir - oo) = exp ( . (12) 





We observe that for all plotted field strengths the curves converge to this limit. With other 
words one can conclude that for large times we obtain the old result given by Schwingers 
formula for weak fields as well as for strong fields. The absolute value of the distribution 
function at r ^ oo is larger for strong fields than it is for weak fields. 

We obtain a similar result for fermion pair creation, see Fig. 3. Due to the different 
amplitudes for fermion and boson pair production, Eqs. (0) and (^), the shapes of the 
curves are slightly different, in particular around zero kinetic momenta. We can compare 
the distribution function for fermions and bosons for a given field strength. Fig. 2 and Fig. 
3. The onset of particle creation (first maximum) is earlier for fermions than it is for bosons. 
Both curves reach the same limit for large times and oscillate with the same frequency for a 
given Eq. The results which we obtain for the boson case are in agreement with those of ||23 



for the case of weak fields, wherein the Markovian limit was employed. But the numerical 
solution introduced in this section allows in addition to consider pair creation for strong 
fields. 
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C. The low density and the Markovian Umit 



1. A closed kinetic equation 

In the previous subsection we have demonstrated how to solve the kinetic equation in 
its non-Markovian form and explored the numerical solutions. In this Section we want to 
discuss the solution of the kinetic equation (|l]) for weak fields where it is possible to apply 
approximations. These are the Markovian approximation and the low density limit . 



The Markovian limit of the non-Markovian equation Eq. (|^) is defined by the neglect of 
memory effects in the source term. Technically this means that we approximate the distri- 
bution function occuring under the integration by /(r') f {''')■ We obtain the following 
kinetic equation 

=[l± 2/f (r)]Si(r) = 5f (r), (13) 

where 5° (r) is the source term in low density limit, Eq. (^. 

Taking into account the initial condition lim /+^(t) = 0, we obtain the following solu- 

T—^ — OO 

tion of the kinetic equation 



/fW=^(l-exp 



±2 r dT'Slir') 



(14) 



This result is exact in the Markovian limit and holds for any time-dependent homogeneous 
electric field. In the lowest order of the expansion of Eq. ([l^ we obtain 



/°(r)=/ dr'Sl{r'). (15) 



This solution is equivalent to the low density limit where effects due to the symmetry 
character of the created particles are neglected in the source term. Note that both the low 
density limit and the Markovian approximation are restricted to weak fields, Eq < 1. 

Before we discuss the numerical results for the Markovian approximation and the low 
density limit, we explore general properties of the distribution function defined in Eq. (|14|). 
One requirement of a kinetic theory is that the distribution function for bosons and for 
fermions must be positive definite for all times and momenta. In order to prove the validity 
of this important property we rewrite Eq. (|T^) as 

fl{r) = I J dr'giir') J dr" gi{r") + \ j dr'gl{r') J dr" gl{r") , (16) 

— oo — oo — oo — oo 

where the functions glf'ir) are given as 

9iV) = w±wl!S®L1j|. (17) 



sin[2e( 



r 



It is easy to transform the integral (ITO) to the following form 
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flir) = \{J dr'glir'))' + 1( / rfr'^l(r'))' . (18) 



From this quadratic representation of Eq. ([Tq ) we can conclude that the distribution function 
is positive definite as it should be, 



flir) > . (19) 



On the other hand this result is very useful for numerical calculations since Eq. ([T8|) just 



requires to perform a single integration compared with Eq. (ITH) in which it was necessary 



to solve a double integration, hence it makes the numerical treatment easier. 

At large times, r — > oo, the distribution function is in the asymptotic regime. In low 
density limit it is easy to show the feature that /o(oo) ^ 0. The distribution function evolves 
from zero at r — — cxd to an asymptotic nonzero value. In the absence of back reaction and 
collisions (or any other damping mechanism) we observe an accumulation effect. In the next 
subsection we will elucidate these properties with numerical results. 



2. Numerical results 

Using the results of Eqs. ([T8|) we can explore the solution of Eq. (p^ for the distribution 
functions in the Markovian limit and in the low density approximation. 

In Figs. 4 and 5 we compare the distribution functions of the non-Markovian solutions 
with the low density limit solutions for fermions and bosons, respectively. In the lower panels 
we have chosen a weak field and observe that the results agree with each other. This was 
expected since for weak fields the absolute value of the distribution function is small and 
consequently the statistical factor [l±2/(r)] does not considerably deviate from one. Hence 
the low density limit is a good approxmation for relatively small field strengths, Eq < 1. 

In the upper panels we have chosen large field strengths. Both for fermions. Fig. 4, and 
for bosons. Fig. 5, the low density limit solution shows a different limit for large times. The 
inclusion of the quantum statistical character in the non-Markovian solution leads for fermion 
pair creation to a suppression and for boson pair creation to an enhancement compared to 
the corresponding low density limit. These plots elucidate two important things. On the 
one hand they demonstrate the influence of the different symmetry character of fermions 
and bosons. On the other hand we see that for strong fields the low density limit solution 
provides wrong results. Technically the reason is clear. The statistical factor [1 ± 2/(r)] 
deviates from one and therefore it has to be included into the kinetic equation. Physically 
this means that the non-Markovian character becomes important. 

We have to distinguish different time scales: the memory time and the production time 
T9| , |22| -p4| . While the memory time has quantum mechanical origin and can be considered 



as the time needed to tunnel the barrier, the production time is the time interval between 
two creation processes. The memory time is now for strong fields of the same order of 
magnitude as the production time. A separation of the time scales, which is necessary 
for the approximation discussed in the previous subsection, is no longer possible. The pre- 
history affects the evolution of the distribution function and therefore memory effects become 
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important. This is taken into account by the time integration over the statistical factor in 
the full non-Markovian equation. 

In the previous subsection we have also discussed the Markovian approximation. It is 
characterized by the neglect of memory, but the different statistical character of fermions 
and bosons is still taken into account. In Fig. 6 we compare the three solutions of the kinetic 
equation for strong fields for bosons (lower panel) and for fermions (upper panel). Besides 
the already discussed non-Markovian solution and the low density approximation, we also 
show the Markovian limit. Needless to point out that for weak fields the Markovian solution 
also agrees with the non-Markovian solution. Although the Markovian solution is in better 
agreement with the correct solution than the low density limit, the error is visible and the 
Markovian limit fails for strong fields. For weak fields the low density limit provides exact 
results since non-Markovian effects disappear. For strong fields it is unavoidable to solve the 
non-Markovian equation. Just in a small band of field stengths of about one the Markovian 
limit is a sensible approximation, better than the low density limit but non-Markovian effects 
are still very small. Hence for small fields, £'0 < 1, we suggest to use the low density hmit, 
for strong fields, £'0 > 1, the non-Markovian solution. 

III. SUMMARY 

We have analytically and numerically explored the solution of a quantum kinetic equation 
describing particle production of boson and fermion pairs. The source term providing the 
creation of pairs is characterized by its non-Markovian character. The time evolution of the 
distribution function depends on the entire pre-history of the evolution. We have numerically 
solved the kinetic equation in its non-Markovian form for both weak and strong fields. We 
compared these solutions with both the low density limit and the Markovian approximation 
for which we found an analytic solution in closed form. 

We observe that for strong fields the results depend on the statistical character of the 
produced particles. We obtain a suppression of the production rate for the creation of 
fermion pairs and an enhancement for boson pairs when we compare the full solution with 
the approximation where memory effects arc neglected. In the domain of weak fields the 
exact solution of the Markovian limit can be approximated by its low density limit being the 
lowest order term of an expansion of the full Markovian solution. As a numerical result we 
obtain that indeed for weak fields the use of the Markovian limit as well as of the low density 
limit is justified. The curves calculated in the Markovian limit, the low density limit and 
the full solution are in good agreement. This means that non-Markovian effects disappear 
for weak fields. Therefore it is sufficient to use for weak fields just the low density limit. 

Furthermore we find that for strong field strengths the Markovian approximation as well 
as the low density limit break down since memory effects become important. For strong fields 
the distribution function has a large value and therefore the statistical factor [1 ± 2/(t)] 
deviates from one. Hence it is transparent that the non-Markovian character appears for 
strong fields. This observation indicates that it is necessary to include memory in order to 
obtain e.g. the correct limit for large times. 

In conclusion we can summarize that for weak fields it is appropriate to use the low den- 
sity limit while for strong fields it is necessary to solve the non-Markovian kinetic equation. 
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In order to give a more complete picture of the physics beyond the Markovian hmit it is 
necessary to study in a next possible step the momentum dependence of the distribution 
functions. 

Although the numerical analysis was performed for a constant field our results quali- 
tatively hold for any time-dependent field. This will be important when a more realistic 
scenario is addressed where the electric field is determined self-consistently. In order to 
incorporate back reactions it is necessary to solve the Maxwell equation that determines the 
electric field via the conduction current and the polarization current due to the creation of 
the charged particle pairs |]TB|,E3[. Furthermore, it would be of great interest to extend this 



approach to the QCD case in order to explore the consequences of the new source term for 
the pre-equilibrium physics in ultrarelativistic heavy-ion collisions. 
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FIGURES 



FIG. 1. The dependence of the source term in low density hmit on the parallel momentum for 
fermion pair creation for a strong field (upper panel: Eq = 1.5) and a weak field (lower panel: 
Eo = 0.5). 

FIG. 2. The boson distribution function as a solution of the non-Mar kovian equation as function 
of time for different field strengths at = 0. 

FIG. 3. Same as Fig. 2 for fermions. 

FIG. 4. The time evolution of the distribution functions of bosons (/_|_) compared with the low 
density limit at = for a strong field (upper panel: Eq = 5.0), and a weak field (lower panel: 
Eq = 0.7), is shown. 

FIG. 5. The time evolution of the distribution functions of fermions (/_) compared with the 
low density limit at = for a strong field (upper panel: Eq = 3.0), and a weak field (lower 
panel: Eq = 0.5), is shown. 

FIG. 6. The time evolution of the distribution functions of fermions (upper panel: Eq = 3.0) 
and bosons (lower panel: Eq = 5.0) within different approximations for strong fields. Only the full 
non-Mar kovian solution provides the correct limit for large times. 
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FIG. 3 



14 




FIG. 4 



15 




FIG. 5 
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